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Clinical Article 



Magnetoencephalography Interictal Spike Clustering 
in Relation with Surgical Outcome of Cortical Dysplasia 
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Objective : The aim of this study was to devise an objective clustering method for magnetoencephalography (MEG) interictal spike sources, and to 
identify the prognostic value of the new clustering method in adult epilepsy patients with cortical dysplasia (CD). 

Methods : We retrospectively analyzed 25 adult patients with histologically proven CD, who underwent MEG examination and surgical resection for 
intractable epilepsy. The mean postoperative follow-up period was 3.1 years. A hierarchical clustering method was adopted for MEG interictal spike 
source clustering. Clustered sources were then tested for their prognostic value toward surgical outcome. 

Results : Postoperative seizure outcome was Engel class I in 6 (24%), class II in 3 (12%), class III in 12 (48%), and class IV in 4 (16%) patients. 
With respect to MEG spike clustering, 1 2 of 25 (48%) patients showed 1 cluster, 2 (8%) showed 2 or more clusters within the same lobe, 1 0 (40%) 
showed 2 or more clusters in a different lobe, and 1 (4%) patient had only scattered spikes with no clustering. Patients who showed focal clustering 
achieved better surgical outcome than distributed cases (p=0.01 7). 

Conclusion : This is the first study that introduces an objective method to classify the distribution of MEG interictal spike sources. By using a hierar- 
chical clustering method, we found that the presence of focal clustered spikes predicts a better postoperative outcome in epilepsy patients with CD. 

Key Words : Epilepsy • Cortical dysplasia • Magnetoencephalography ■ Hierarchical clustering ■ Surgical outcome. 



INTRODUCTION 

Cortical dysplasia (CD) is the second most common patho- 
logic entity in surgically treated epilepsy 4 ' 8,10 ' 11 - 13,171 . However, 
the postoperative seizure outcome is relatively limited when 
compared outcome with epilepsy associated with hippocampal 
sclerosis or tumor 4,718 '. Recently, magnetoencephalography 
(MEG) has emerged as a useful tool for localization of the epi- 
leptogenic zone. MEG studies, which have attempted to deter- 
mine the relationship between interictal spike source patterns 
and postoperative seizure outcome, have mostly dealt with 
children and their patient groups were mixed with several dif- 
ferent pathologies 9,15,251 . Moreover, those studies used a subjec- 
tive clustering method to determine clusters; none of the stud- 
ies used a statistically based clustering method. One MEG 



study recruited children who had a single pathology of CD 251 . 
However, because they were interested in influences of com- 
plete or partial resection of areas containing MEG spike sourc- 
es on surgical outcome, they did not report an association of 
clustering pattern with surgical outcome. There is one MEG 
study with adult epilepsy patients associated with CD, but that 
study included only four patients 141 . Therefore, based on the 
above studies alone, it is difficult to determine which features of 
the MEG spike source have prognostic value in CD-associated 
epilepsy surgery. 

In this study we developed an objective clustering method to 
classify MEG spike source patterns, and we determined the 
prognostic value of that new clustering method in adult epilep- 
sy patients with CD. To that end, we undertook a retrospective 
analysis of MEG data with respect to postsurgical outcome. 



• Received : May 23, 201 2 • Revised : July 25, 201 2 • Accepted : November 22, 201 2 

• Address for reprints : Chun Kee Chung, M.D., Ph.D. 

Department of Neurosurgery, Seoul National University Hospital, 1 01 Daehak-ro, Jongno-gu, Seoul 1 1 0-744, Korea 
Tel : +82-2-2072-2352, Fax : +82-2-2072-0806, E-mail : chungc@snu.ac.kr 

• This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.0rg/licenses/by-nc/3.0) 
which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited. 



466 



MEG Interictal Spike Clustering W Jeong, et al. 



MATERIALS AND METHODS 

Patient population 

Initially included were 58 patients with intractable epilepsy 
with histologically proven CD who underwent MEG examina- 
tion and surgical operation between 2005 and 2010 at our hos- 
pital. Of that total, 15 patients with less than a 1 year follow-up, 
one patient younger than age of 15 years at surgery and 8 pa- 
tients who had fewer than 5 MEG interictal spikes were exclud- 
ed. Also excluded were 8 patients with medial temporal lobe 
epilepsy and one patient who had atrophy. As a result, 25 pa- 
tients (mean age at surgery=29.3 years; SD=8.3; 11 females) 
were included in our retrospective study. The mean postopera- 
tive follow-up period was 3.1 years (SD=1.2 years). All patients 
provided informed consent. This study was approved by the In- 
stitutional Review Board of Seoul National University Hospital. 

Histological classification 

We classified focal CD (FCD) according to Blumcke's crite- 
ria 21 . FCD type I is characterized by radial (FCD type la) or tan- 
gential (FCD type lb) dyslamination of the neocortex. FCD 
type II is characterized by cortical dyslamination and dysmor- 
phic neurons without (type Ha) or with (type lib) balloon cells. 
FCD type III was not included because, by definition, this type 
is associated with other pathologies. 



cardiogram were simultaneously recorded in order to monitor 
eye movement and cardiac artifacts. The location of the subject's 
head in the MEG sensor space was determined by measuring the 
magnetic signals from four head position indicator coils placed 
on the scalp. The coil locations with respect to the three anatomi- 
cal landmarks, the nasion and two preauricular points, were 
identified by a three dimensional digitizer (FASTRAK™, Polhe- 
mus, Colchester, VT, USA). The obtained information was used 
to co-register the MRI, MEG, and head coordinates so that the 
sources of the MEG signals could be superim-posed on the sub- 
jects' MR image. After the recording, a signal-space-separation 
method was applied in order to reduce environmental and bio- 
logical noises 21,22 '. Each MEG interictal spike was visually inspect- 
ed and localized using the single equivalent current dipole (ECD) 
method in Neuromag software (Elekta Neuromag, Helsinki, Fin- 
land). The ECDs were identified by a least-squares search in a 
spherical volume conductor model by using subsets of 20-50 
channels surrounding the maximum signals. Only ECDs with a 
goodness-of-fit value >85% and confidence volume <3 mm 3 
were accepted. Accepted ECDs then were superimposed on the 
patients' individual MR image. 

MEG spike clustering 

Hierarchical clustering analysis was adopted for MEG spike 



Magnetic resonance imaging 

All patients were examined using ei- 
ther a GE 1.5 or 3 T MRI unit (GE Ho- 
rizon Echospeed) or a Siemens 1.5 T 
scanner (Siemens Avento system, Erlan- 
gen, Germany) both before and after 
surgery. The images were reconstructed 
offline and transferred to a UNIX work- 
station (HP C3750). 

MEG acquisition & spike source 
localization 

Cortical magnetic activity and electro- 
encephalographic (EEG) activities were 
recorded simultaneously using a whole- 
head MEG system (VectorView™, Elek- 
ta Neuromag Oy, Helsinki, Finland) 
consisting of 306 channels arranged in 
triplets of two planar gradiometers and 
one magnetometer. Spontaneous MEG/ 
EEG signals were acquired about 60 min 
(four 15 min epochs of data) while the 
patients were lying in a supine position 
with eyes closed in a magnetically shield- 
ed room. The sampling rate was 600 Hz 
with a bandpass filter of 0.1-200 Hz. A 
bipolar electrooculogram and an electro- 
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Fig. 1 . Hierarchical dendrogram (A) and equivalent current dipoles on MR images (B) for Patient #1 0. 
The individual spike sources (n=20) are represented along with the vertical axis. The ^-shaped lines 
connect the spike objects that are combined at each aggregation stage. 
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source clustering using SPSS 19.0 software (SPSS, Chicago, IL, 
USA). Euclidean distances between ECDs in each patient were 
calculated as input values. The hierarchical grouping starts by 
defining a number of clusters equal to the total number of 
ECDs, N. In the first grouping step the number of clusters is re- 
duced to N-l by combining the two closest ECDs using a com- 
plete linkage method. In successive stages the newly formed 
groups are linked to other ECDs until the entire ECD popula- 
tion is finally linked together into one cluster. To make an ap- 
propriate number of cluster groups, a dendrogram tree, which 
shows the result of the hierarchical clustering analysis, was cut 
by using a two-step criteria. First, if the ECDs projected onto 
the MR image were shown as clearly different clustered loca- 
tions, we divided the dendrogram where the same number of 
clusters are grouped (Fig. 1). Second, we calculated the stan- 
dard deviation (SD) of the Euclidian distances of ECDs within 
the same cluster. The dendrogram tree was then cut where any 
pair of ECDs distances within the same cluster did not exceed 
the 2.5 SD. The 1 SD should not exceed 10 mm. More than 3 
spikes in one group was considered a cluster. 

Surgical procedures and postsurgical seizure outcome 

Extraoperative intracranial recording was done for all 25 pa- 
tients to delineate the ictal-onset zone and to obtain functional 

Table 1 . Summary of patient profiles 



mapping. The extent of invasive intracranial electrode implant 
was determined as follows 4) . If MRI showed a visible lesion, in- 
vasive intracranial electrodes were placed around the lesion. If 
MRI showed no visible lesion, invasive intracranial electrodes 
were placed to encompass the ictal-onset zone delineated by 
scalp video-EEC monitoring. When noninvasive results were 
discordant, invasive intracranial electrodes were implanted bi- 
laterally as needed. After the invasive study, video-EEC moni- 
toring was done until at least three typical seizures recorded. 
Resection included the whole extent of the visible lesion, if 
present, or the whole ictal-onset zone bound by sulci or func- 
tional cortex. In addition, underlying white matter was resect- 
ed, if feasible. The surgical outcome was classified according to 
modified Engel criteria 5 '. 

Statistical analysis 

All analysis was performed by using SPSS 19.0 software (SPSS). 
We used the X 2 test or, where appropriate, Fishers exact test. A 
value ofp<0.05 was considered to be statistically significant. 

RESULTS 

Postoperative seizure outcome was Engel class I in 6 (24%) 
patients, class II in 3 (12%), class III in 12 (48%), and class IV in 
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4 (16%) patients. Table 1 presents a summary of the patients' 
profiles. 

MEG 

Twelve of 25 (48%) patients showed 1 cluster, 8 (32%) showed 
2 clusters, 3 (12%) showed 3 clusters, 1 (4%) showed 5 clusters, 
and 1 (4%) had scattered spikes with no clustering. Of the pa- 
tients who showed 2 or more clusters, 2 patients' clusters were 
located within the same lobe and 10 patients' clusters were in 
different lobes. We defined the presence of 1 cluster and >2 
clusters within the same lobe as focal clustering and >2 cluster 
in different lobes and scattered spikes as distributed clustering. 
Patients with focal clustering had a better surgical outcome 
than those with distributed clustering (p=0.017, Engel I vs. II- 
IV) (Table 2). 

We visually checked whether the MEG clustered area was re- 
sected regardless of its extent (complete or partial) by project- 
ing the ECDs on postoperative MR images. The area that con- 
tained a MEG cluster was removed in 5 of 6 (83%) cases in the 
Engel class I group, 2 of 3 (67%) in Engel II, 5 of 12 (42%) in 
Engel III, and 1 of 4 (25%) in the Engel IV group. Patients who 
had resection of a MEG cluster tended to show a better surgical 
outcome (p=0.097) (Table 1). 

Histological subtypes and resection location 

Of the 25 patients, 19 (76%) had type I and 6 (24%) had type 
II CD (Table 2). With respect to postoperative seizure outcome, 
there were no statistically significant differences were found be- 
tween patients having type I CD and type II CD (p=0.074). Of 
the 25 patients, 14 (56%) underwent 
frontal lobe resection, 8 (32%) tempo- 
ral lobe resection, and 3 (12%) occipital 
lobe resection (Table 2). There was no 
postoperative seizure outcome differ- 
ence among patients with different lobe 
resection areas (p=0.278). 



adopted to classify MEG spike sources in 25 adult patients with 
histologically proven CD, who underwent surgical resection for 
intractable epilepsy and had at least 1 year of follow-up. The re- 
sults of the MEG spike source clustering method were then 
tested for their prognostic value toward surgical outcome. The 
14 (56%) patients in the focal clustering group achieved better 
outcome than the 1 1 (44%) in the distributed clustering group. 
No statistical differences were found according to the CD type 
and MRI findings with respect to surgical outcome. 

Hierarchical clustering method 

Some researchers have aimed at determining the relationship 
between MEG interictal spike source patterns and postopera- 
tive seizure outcome in their studies 9,15,251 . However, the fact that 
each researcher used different criteria in classification of spike 
sources to describe their results makes it difficult to evaluate the 
prognostic values of MEG in epilepsy. Moreover, all of these 
clustering methods are somewhat subjective in that none used a 
statistically based clustering method. The hierarchical clustering 
method which we adopted in this study is a widely used data 
analysis tool in many fields, including machine learning, data 
mining, bioinformatics, and image analysis. It is a simple and 
objective method compared to other clustering methods used 
by previous researchers. To adopt this method in MEG spike 
source clustering, we simply entered the Euclidean distance 
matrix between ECDs as an input value. Subsequently, the ob- 
tained dendrogram depicts the clustering results by linking the 
closest ECDs together. By using this clustering method, we 
could clearly separate sub-groups of either focal or distributed 



Table 2. Surgical outcome with respect to MEG cluster, CD type, resected lobe, and MR findings 



Engel class 



III 



IV 



MRI findings 

Nine (36%) of 25 patients showed a 
visible lesion of focal CD on their MR 
images (Table 2). No statistical differ- 
ences were found between normal and 
abnormal MRI findings with respect to 
surgical outcome (p=0.363). 

DISCUSSION 

In this study, we developed an objec- 
tive method of MEG interictal spike 
source classification, and identified the 
prognostic values of MEG dipole clus- 
ters in adult epilepsy patients with CD. 
A hierarchical clustering method was 
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cluster according to the dendrograms clustering patterns. If other 
researchers who are interested in MEG spike source clustering 
use this hierarchical clustering method, we can expect an in- 
crease in the comparability of results among studies in the future. 

MEG cluster with respect to surgical outcome 

Studies focused on MEG clustering with respect to surgical 
outcome have suggested that focal clustering of MEG spike 
sources and resection of that area may predict a favorable postop- 
erative outcome. However, most of these studies have dealt with 
children or with mixed groups with various pathologies 6 ' 9 ' 15 ' 16,25 '. 
Although there is one study dealing with adult epilepsy patients 
associated with CD, only 4 patients were included in their study 14 '. 
We studied an adult group of 25 patients with histologically 
proven CD. This is the largest MEG epilepsy study with the sin- 
gle pathology of CD in an adult group with more than a 1 year 
follow-up period. 

As a result of the MEG spike source clustering method, 14 
(56%) of our patients showed focal clustering and 1 1 (44%) 
showed distributed clustering. There are statistically significant 
differences in postoperative seizure outcome between these two 
groups : the focal clustering group achieved a better outcome 
than the distributed clustering group. In addition to clustering 
type, we also compared the clustered area with the surgically re- 
sected area as seen on postoperative MR images. The area that 
contained the MEG cluster was removed in 5 of 6 (83%) patients 
in the Engel I group, 2 of 3 (67%) in Engel II, 5 of 12 (42%) in 
Engel III, and 1 of 4 (25%) in the Engel IV group. Four of the 5 
(80%) patients in Engel III had another cluster (other than the 
resected cluster) in a location remote from the actual resection 
area compared with only 1 of 7 (14%) patients with another 
cluster in the Engel I and II groups. One Engel IV patient (#25) 
is the only patient in our group who underwent multiple area 
resection. Additionally, he had an awake-lesionectomy with in- 
tra-operative sensory, motor, and language mapping to avoid 
injury to the eloquent cortex. Multiple area resection and hin- 
drance of complete resection caused by the eloquent area were 
considered as a bad prognostic factor 12 '. Failure to achieve a sei- 
zure free outcome despite resection of the focal clustered area in 
Patient #25 may result from these reasons. 

On the basis of our results, we propose that the presence of 
focal clustered spike sources and resection of that area predict a 
favorable surgical outcome, whereas the presence of distributed 
clustering predicts a poor surgical outcome, in adult epilepsy 
patients with CD. 

MEG cluster with non-lesional cases 

To draw firm conclusions about the relationship between 
MEG clustering and surgical outcome, explanation of the pa- 
tients who had poor outcome despite the presence of focal clus- 
tering is necessary. This is especially important for the patients 
who showed no abnormality on their MR images. Many studies 
agree on the prognostic value of MRI -visible lesions 1,19,23,26 ', but 



not many studies have elucidated the prognostic factors of the 
patients with normal MRI findings, especially in relation with 
associated electrophysiological characteristics 16 '. Of our patients 
who showed no abnormal features on their MR images, 3 pa- 
tients showed a poor surgical outcome (Engel class IV) despite 
the presence of focal MEG spike clustering. One patient (#25) 
had an awake-lesionectomy with multiple resection of frontal 
lobe as stated previously, while two other patients (#22, #23) had 

I cluster located in the posterior part of the resection area. We 
suggest that if the area containing MEG spike sources were re- 
sected, our non-lesional, single clustered patients might achieve 
favorable seizure outcome. Further intensive study with non-le- 
sional patients is needed to confirm this suggestion. 

Histological subtypes and MR findings 

With regard to postoperative seizure outcome and CD type, 
no statistically significant differences were found. Even though 
there is a debate about surgical outcome depending on CD type, 
recent studies have concluded that type II CD patients have bet- 
ter outcome than type I patients 3,4,20,24 '. The discrepancy be- 
tween previous studies and our results may due to the limited 
number of patients. Tendency to show better outcome in type 

II CD in our results support this idea (p=0.074). Even though 
many studies report better outcome in lesional cases, no surgi- 
cal outcome differences were observed between patients with 
normal and abnormal MRI findings in our patient cohort. This 
could be explained by the fact that the small number of patients 
having abnormality on MR images in our study makes it diffi- 
cult to reveal statistically significant results. Only 9 of 25 pa- 
tients showed lesions on preoperative MR images. Furthermore, 
5 of 6 who achieved Engel class II-IV had divergent presurgical 
evaluation results, including PET, ictal SPECT, video-EEG mon- 
itoring, invasive EEG, MEG, and MRI. Previous studies have re- 
ported that divergent presurgical evaluation results predict 
poor outcome 13,16 '. Thus, the limited number of patients and the 
presence of divergent presurgical evaluation results may play 
roles in the lack of improved outcome in the lesional cases in 
our study. 

Limitation 

Its retrospective design in a single institute is a major caveat 
which should be taken into account in exploiting the results of 
the present study. Nevertheless, in terms of the prognostic val- 
ues of the MEG interictal spike cluster, this is the largest MEG 
study that has dealt with adult epilepsy group who had the sin- 
gle pathology of cortical dysplasia. In addition, this is the first 
study which introduces an objective method to classify the dis- 
tribution of MEG interictal spike sources. 

CONCLUSION 

In this study, a hierarchical clustering method was adopted 
for MEG interictal spike source clustering in adult epilepsy pa- 
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tients with CD. By using this method of hierarchical clustering, 
we were able to delineate the value of using clustering to predict 
surgical outcome more clearly. The presence of focal clustered 
MEG interictal spike sources and subsequent resection of the 
clustered area predicted a favorable postoperative outcome. In 
contrast, the presence of a distributed cluster predicted a poor 
surgical outcome. 
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